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Recently, the Gravitational-Search Algorithm (GSA) has been presented as a 
promising physics-inspired stochastic global optimization technique. It takes 
its derivation and features from laws of gravitation. This paper applies the 
GSA to design optimal controllers of a nonlinear system consisting of a 
doubly-fed induction generator (DFIG) driven by a wind turbine. Both the 
active and the reactive power are controlled and processed through a back-to- 
back converter. The active power control loop consists of two cascaded 
proportional integral (PI) controllers. Another PI controller is used to set the 
q-component of the rotor voltage by compensating the generated reactive 
power. The GSA is used to simultaneously tune the parameters of the three 
PI controllers. A time-weighted absolute error (ITAE) is used in the objective 
function to stabilize the system and increase its damping when subjected to 


different disturbances. Simulation results will demonstrate that the optimal 
GS A-based coordinated controllers can efficiently damp system oscillations 
under severe disturbances. Moreover, simulation results will show that the 
designed optimal controllers obtained using the GSA perform better than the 
optimal controllers obtained using two commonly used global optimization 
techniques, which are the Genetic Algorithm (GA) and Particle Swarm 
Optimization (PSO). 
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1. INTRODUCTION 

Generated electric power from wind is growing at an increasing annual rate worldwide in the last 
decade as its cost has decreased considerably. Widespread use of wind power can largely reduce production 
costs, reduce carbon emissions, contribute to resolving global warming calls, and reduce the dependency on 
fossil fuels [1]. 

With the recent developments in power electronics technology, DFIG driven by a wind turbine [2] 
are becoming very popular and their worldwide market share is increasing quickly. DFIG has been widely 
used in grid-connected variable-speed constant-frequency wind energy conversion systems (WECS) due to 
their advantageous features such as higher efficiency, lower cost, wide-range operations, capability of 
producing reactive power, better energy capture, smoother operation, and flexible control. Therefore, with the 
increasing utilization of wind power in the interconnected power network, improving DFIG transient 
performance is becoming more significant in power system operations. 

Tuning the parameters of a classical PI controller of a grid-connected DFIG-based WECS is 
complicated due to the high nonlinearity and complexity of the system. Optimization techniques are applied 
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to determine the optimal parameters of a classical PI controller to enhance the performance of the WECS as 
its response change certainly with time. In [3] an adaptive fuzzy gain scheduling of PI controller is presented. 
The PI controller parameters are optimized online using an adaptive mechanism based on fuzzy logic. In [4] a 
non-dominated sorting GA alongside discrete Fourier transform is implemented to quantify the damping 
characteristics of the rotor-side current due to system disturbance. In recent years, different evolution 
algorithms [5], [6] including mean variance optimization (MVO) algorithm [7], GA [8], PSO [9-12], bacteria 
foraging optimization (PFO) algorithm [13], hybrid artificial neural network [14] and Gbest guided artificial 
bee colony algorithm [15], [16] are applied to optimize the parameters of the PI controller of the rotor side 
converter (RSC) and hence improve the damping of oscillatory modes in DFIG-based WECS. In [17] a PSO 
is presented to generate an on-off controller to ensure robustness and quality of the energy produced, wherein 
[18] a nonlinear rotor-side controller is designed based on H, optimal control theory to demonstrate the 
synthesis of a maximum power point tracking (MPPT) algorithm. Authors in [19] presented an optimal 
tuning of the MPPT algorithm and the active disturbances rejection controller to reduce the impact of the 
change on the characteristics of DFIG-based WECS. In [20], authors presented a computational intelligence 
control method to optimize the performance of the back-to-back converter and to regulate the low voltage 
ride-through capability of DFIG-based WECS. In [21], the DFIG-based WECS small signal stability was 
analyzed and their impacts on the eigenvalues of a DFIG in a single machine infinite bus system was 
investigated. In [22], the states of DFIG are used to obtain the optimal feedback control with utilizing the 
rotor currents beside the rotor voltages. Fuzzy logic nonlinear controller and particle swarm optimization are 
used in [23] to handle the nonlinearities in DFIG and to optimize the controllers’ gains. In [24], the authors 
optimized the performance of the grid side inverter (GSI) and RSC of grid-connected DFIG to minimize the 
tracking error based on discrete-time PSO inverse optimal controller. Differential geometry theory is used in 
[25] to control linearized DFIG-wind turbine model based on linear quadratic regulator. To enhance the 
DFIG during and after power system transients by optimizing the PI-controllers of the RSC, PSO is used in 
[26] and GA in [27]. In [28], the transients of the DFIG are optimized by minimizing the sum squared error 
deviation of the dc link voltage and the generator speed deviation. 

GSA is one of the recent heuristic algorithms, which has been introduced as a physics-inspired 
method in 2009 [29], [30]. It takes derivation and the main characteristics of heuristics from “Newton’s law 
of gravitation”; “Every particle in the universe attracts every other particle with a force that is directly 
proportional to the product of their masses and inversely proportional to the square of the distance between 
them” [31]. The GSA updates gravitational and inertia masses with the help of heavy masses and finds the 
optimum [29-31]. In the last few years, GSA has been successfully applied to solve world engineering 
optimization problems and demonstrate high quality performance in solving these problems [32-35]. 

In this paper, GSA is used to design optimal controllers for improving transient response of a 
nonlinear DFIG-based WECS. The optimal controllers are tested on different common types of disturbances. 
The results are compared to the results obtained using the two well-recognized global optimization 
algorithms, which are GA and PSO, to show the effectiveness of using GSA to attain a global optimal 
solution of the design problem under study. 


2. MODEL OF THE DFIG-BASED WECS 

Figure 1 depicts a grid-connected DFIG-based WECS. The system shown includes a wind turbine, 
an induction generator, a back-to-back converter (RSC and GSI) [36], and the utility-grid. The ac/dc and the 
dc/ac conversions are used to control the bidirectional power delivered from/to the rotor circuit and to/from 
the grid, respectively. The dynamic electrical and mechanical equations are summarized below. 

The induction generator considered in this paper is a wound rotor induction machine. If the 
magnetic branch is assumed linear, the dynamic modeling of DFIG in arbitrary reference frame can be 
written as seen from the stator side as follows: 


vgs = "sigs + @Ads + PAgs (1) 
Yas = "sas 7 OAgs + Pras (2) 
Vr =T'r Ugrt(o— Op )A' qg PA ar (3) 
ae =r'y i' g Tr JA'qr +PA' dr (4) 
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Ags = Usigs +M Cigs +i'gr ) 
qs 1s'qs qs tar (5) 


GEAR Box 


GENERATOR SIDE GRID SIDE 


Figure 1. General configuration of grid-connected DFIG-based WECS 


Aas = sigs +M lgs + ay) 6) 
Agr = Ltr ï'gr+M Cigs +i'gr) (7) 
Bar = Er egy tM Uys +8 Gy) (8) 


where y > Ags and 4 ds ae the stator voltages, currents, and flux linkages in qd-reference 


v i i 
9S? d? “98? ds 


frame, respectively. Vi vdr» iqr>idr> Agr 


reference frame, respectively. œ is the arbitrary angular speed, w, is the electrical angular speed of the 


and Adr are the rotor voltages, currents, and flux linkages in qd 


rotor, , is the derivative operator, rọ is the stator windings resistance, r', is the rotor windings resistance, 


M =3/2Lns.» Ls is the stator magnetizing inductance, Iris the rotor leakage inductance, and Lj, is the 


stator leakage inductance. 
If P indicates number of poles, the electromechanical developed torqueTy can be expressed as: 


3P, , 
Ta = 55 Gasiqs — Aqsias) (0) 


The dynamics of the mechanical shaft can be expressed based on the torque applied on the shaft Ty 
from both the generator end and the turbine end as follow: 


Tsh =Tr +Jr 


dt dt (10) 


where T, is the applied mechanical torque at the turbine side, T, is the electrical torque seen at the 
generator side, J; and Jm are the inertia of the turbine and the generator, respectively, and œz and @y, are 
the angular speed of the turbine and the generator, respectively. The existence of the speed difference on the 
shaft ends and the compliance coefficient of the shaft Ksp can determine the shaft-torsion (Ag — 0z ) as 


in (11): 


Tsh 
bo —O, =((@n-0,)-dt == 
gle Nene a a 
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As shown in Figure 2, both the active and the reactive power are controlled through the back-to- 
back converter. The active power control loop, which sets the d-component of the rotor voltage y a consists 
r 


of two cascaded PI-controllers; the first one is used to regulate the rotor speed and to set the reference active 


power to the second PI-controller which is used to compensate the generated active power p, . Another PI 


controller is used to set the q-component gr by compensating the generated reactive power g,. The 


mathematical description of the PI controllers is given by: 


AC GRID GENERATOR SIDE 


' A 
V dr V gr 
SPEED CONTROL ACTIVE POWER CONTROL 
> K pps) 
Kips >| Sx(adt 
y(t) “ x(t) i(Ps) 
pt REACTIVE POWER CONTROL 

>| Kros) 

Qs ref + gi Ly + 

2 wO Kios) w(t)dt 2 


o| 
Figure 2. Block diagram of the grid-connected DFIG-based WECS with its associated control system 


; t 
Vdr (t) = K pP x(t) + KiP, Jx(ojdo 
0 (12) 


t 
YO =K pw g(t) + Kio | Og (odo 
o (13) 


r t 
var ©) = K pos WO + Kio, f wodo 
o (14) 


The constants K pœ- Kio» K PPs ; Kip, ; K bOs ; Kios are the proportional gains and integral 


gains of the PI controller as shown in Figure 2. Each gain of them is responsible for doing a specific task 


made on the input deviation signal; proportional gains ( K po’ Kp Ps > K po is ) are concerned with system 


stability; the integral gains ( Kjq. Kip,» KiQg ) are used of the steady-state error. In order to restrict the 


level of output of the PI controlled speed fluctuation during transient conditions, limits are imposed as shown 
in Figure 2. 

The parameters of the control system to be optimally tuned are the six PI controller constants 
( Kpo Kio» K PPs $ KpQs ; KiQs ). In control system analysis and design, performance indices are 


used to assess a system behavior. This behavior is judged to be optimal when the performance induces reach 
minimum values and all system constraints are satisfied. In this paper, the integral of the time-weighted 
absolute error (ITAE) measure is used, which is given as in (15). 
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ts 
ITAE = Ĵt | e(t) | dt 


0 (15) 


where ft, is the simulation time, and e(t) is the error signal which represent the difference between a 
reference value and a measured variable. The error signal is selected as the speed error signal 
Celt) = @m,ref — 2m). 

The optimal values of controller gains are optimally tuned solving the following generalized 
optimization problem: 


Minimize ITAE (K po Kia» K pQg KiQs K pP, > KiPs) (16) 
Subject to 


pgr < Kpo < Kmax 
po po 


K™in < Kjo < K™max 
i@ io 


min <K < gmax 
PQS PQs PQS 


KMn < Kio, < KMX 
iQ g iQs iQs 
gmin < KpP; < K max 
pPs PPs 


min < g.p < pmax 
Kip, $ Buby $ Bip. 


3. GRAVITATIONAL SEARCH ALGORITHM 

GSA is a physics-inspired optimization method as it is derived and takes derivation and the main 
features of heuristics from Newton's laws of gravity. It depends on the gravity and mass interactions [29-36]. 

Based on Newtonian theory, three masses are assigned to a body, which describe its quite different 
properties; inertial mass(My7), passive gravitational mass (Mp), and active gravitational mass 
(M 4) [29-31]. 

The mass Mņ is the quantity presented in Newton's second law, which is a measure of a body's 
inertia. It is a measure of body's resistance to change in motion when a force is applied on the body, i.e., the 
body opposes that change through its inertial mass. The larger M7, the harder it is to accelerate for a given 
strength of force. Mp is a measure of body's reaction to a gravitational field, i.e., it measures a body's 
response to being placed in a gravitational field, which implies that if M p is placed at some point that has 
gravitational potential, then it will experience a force on it. M Ą is a measure of body's source strength for 


producing a gravitational field by itself [29-32]. 

Newton’s law of gravity indicates that the force between two particles is directly proportional to the 
product of their masses and inversely proportional to the square of the distance between them, which can be 
written as [29]: 


F =G M]M21/R° (17) 
where F is the magnitude of the gravitational force, G is the gravitational constant, Mq and Mp) are the 


masses of the first and second particles, respectively, and r is the distance between them. 
Newton’s second law indicates that a particle's acceleration ( a) depends on the force applied to the 
particle and the inertial mass (M7 ) of the particle, which is expressed as: 


a=F/My] (18) 
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As in the case of physics where performance of objects (particles) is measured by their masses, the 
agents in this method are attracted to each other based on gravitational force, and their tendency is towards 
heavier masses, and the heavy masses move in a manner much slower than the lighter ones, and this 
represents the exploitation behavior. 

The GSA algorithm depends on the fact that every mass (agent) has four specifications: position, 
inertial mass, active gravitational mass, and passive gravitational mass. The algorithm updates gravitational 
and inertia masses with the help of heavy masses and finds the optimum. This artificial world of masses 
obeys Newton’s law of gravity and motion. 

To present the GSA algorithm, consider a system with N agents (masses). The position of the ith 
agent is defined by [29]: 


yrsa za ex’) ;i=1,2,...,N (19) 


l 


d 


where x; is the position of ith agent in the dth dimension. The force acting on mass ; from mass j at a 


specific time z , is defined as follows [29]: 


M Pi ©M aj) 


Fe =G(t) 


f (x4 () xf 0) 


R.. (t) +E 
= (20) 
where M p;(t) is the passive gravitational mass of agent ;, M Aj) is the active gravitational mass of 


agent J, G(t) is the gravitational constant at time z, s is a small constant, and Ri) is the Euclidian 


distance between two agents ; and j [29]: 


Rj =o X ;O| (21) 


The total force that acts on agent ; in a dimension d is a randomly weighted sum of the dth 
components of the forces exerted by other agents, which is expressed as [29]: 


N 
FIO = È rand; FIO 
j=l, j#i (22) 
where rand j is a random number in the interval /0,]] , which is introduced the equation to give a stochastic 
characteristic to the algorithm. If M y; (t) is the inertial mass of ith-agent, then by the law of motion, the 
acceleration in direction d of the agent ; at time +, is given by [29]: 


d d 
af (t)= FS (t)/M 7 (0) (23) 


The next velocity of an agent is considered as a fraction of its current velocity added to its 
acceleration. Therefore, its position and its velocity could be calculated as follows [29]: 


vd @t+D=7 xvid (t) + ad (t) (24) 


dispar? viita 
xf at= av C4) as 
where r; is a uniform random variable in the interval /0, 1] , which is employed in the equation to give a 


randomized characteristic to the search. The gravitational constant G is a function of the initial value ( Go) 


and time (;), i.e., it is initialized at the start and will be decreased as time passes to control the search 
accuracy. This constant is expressed as [29]: 
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Gt) =G(Go,t) =Gp AD (26) 


where g is a constant specified by a user, z is current iteration and 7 is the total number of iterations. 
Fitness evaluation is used in this algorithm to calculate gravitational and inertia masses. A more 
efficient agent has a heavier mass, which indicates that better agents have higher attractions with slower 


walks. The values of gravitational and inertia masses, which are assumed equal, are calculated using the map 
of fitness. The gravitational and inertial masses are updated by the following Equations [29]: 


N 
Mismi O/ VmjO 


j=l (27) 
my (t) = (FAO -WO) BO -W) (28) 
Maj =Mp; =M 7 =M;; i=1,2,...,N (29) 


where F;(t) represents the fitness value of the agent i at time +. For a minimization problem, W(t) and 
B(t) are defined as follows [29]: 


BO)= min FiO) 
Je{l,...N} (30) 


W(t)= max Fi) 
Je{l,...N} 31) 


For a maximization problem, W(t) and B(t) are defined as follows [29]: 


BO= max FiO 
Jefl,...m} (32) 


WO= min FO 
je{l,...m} (33) 


Equation (25) can be modified by reducing the number of agents with lapse of time to perform a 
good compromise between exploration and exploitation, where only a set of agents with bigger mass apply 
their forces to the others. Let Kp is the set of the first x agents with the best fitness values and biggest 


masses, (25) can be written as [29]: 


N 
FL = rand; FIO 
jEKp,j=+i (34) 


The function Kp is a function of time, which has an initial value of Kọ at the beginning and 
decreases as time passes, which implies that all agents apply the force at the beginning, and as time passes 
Kp is decreased linearly. Finally, just one agent with the best fitness values and biggest masses apply force 
to the others. 


4. SIMULATION RESULTS 
This section presents the simulation results of the DFIG-based WECS using the GSA-based 
coordinated optimal tuning of PI controllers. Per-unit system parameters are summarized in Table 1. In order 
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to indicate the effectiveness of the GSA to tune the parameters of PI controllers, the results will be compared 
with the results obtained using GA and PSO as shown Table 2. The GA and PSO are selected in the paper for 
comparisons as they are becoming very common for solving the optimization problems in different 
disciplines, essentially in the design of the controllers because of being robust in finding an optimal solution 
and capable to attain a near-optimal solution close to a global minimum. 


Table 1. Parameters of the Simulated System 


Stator windings resistance Is 0.05 [p.u.] 
Rotor windings resistance ry 0.02 [p.u.] 
Stator magnetizing inductance Lms 5.00 [pu] 


Rotor leakage inductance i 1 0.10 [pu] 


Stator leakage inductance Eis 0.10 [p.u.] 
Compliance coefficient of the shaft K; h 0.40 [p.u./rad] 
Inertia of the turbine Hz; 2.50 [sec] 
Inertia of the generator H g 0.50 [sec] 


Two disturbances are discussed in this paper to show the advantages of using the GSA to tune the PI 
controller over the other algorithms (GA and PSO). Two performance indices are used to compare between 
the results. The first one is the objective function (ITAE) and the second one is the settling time for 5% 
criterion. The comparisons between the three algorithms are listed in Table 2 and Table 3. 


Table 2. Optimal Controller Parameters of GSA, GA, and PSO 


Fault Disturbance Reference Speed Disturbance 
GSA GA PSO GSA GA PSO 
ITAE 0.10585 0.28642 0.61399 1.3265 1.38187 4.63480 
Kyo 20.00 16.881 19.146 20.000 18.954 17.629 


00 35 40 0 94 59 
Kio 20.00 16.384 14.289 20.000 19.036 19.027 
00 74 63 0 54 92 
Kyo 0.659 0.5154 0.4850 10.131 0.0826 0.8635 
10 7 1 7 3 0 
Kio 0.050 3.9751 17.598 20.000 3.7749 6.4840 
00 7 22 0 6 3 
Kpps 0.050 0.2964 0.0500 16.396 9.2915 
00 0.1587 6 43 5 


Kips 0.641 3.1956 14.761 15.920 14.143 18.090 
30 2 53 4 75 68 


Table 3. Calculated Settling time [in seconds] 


Fault Disturbance Reference Speed Disturbance 
GSA GA PSO GSA GA PSO 
Var 0.82 1.19 1.38 2.40 2.71 2.82 
Var 0.10 023 085 156 179 185 
Om 0.04 0.06 0.11 1.04 1.42 1.54 
Ps 0.83 0.96 1.13 2.34 2.6 2.68 


Q, 0.22 0.29 0.94 0.47 0.72 0.72 


By solving the optimization problem using GSA, GA and PSO, the optimal controller parameters of 
the system under study are obtained and summarized in Table 2. As can be seen, the values of ITAE 
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associated with the GSA-based optimal controllers are better than those obtained using the GA-based and 
PSO-based optimal controllers. This fact will be reflected in the following response subsections, Figures and 


Tables. Only five variables (v' gr? Vdp? P;, Qs and @pm ) are presented in this paper while the response of 


the other variables can be estimated from the presented results. 


4.1. Three-phase to Ground Fault 

Three-phase to ground fault was applied after the system reaches to the steady state (t=/0 s) and 
cleared after 30 ms. From Table 2, 60% and 80% reduction of JTAE is achieved using GSA relative to GA 
and PSO, respectively. Figure 3 and Table 3 show that v' dr and v'gr can reach the new steady state much 


faster when GSA is adopted. Figure 4 and Table 3 manifest that using GSA helps the rotor speed to damp in 
40 ms which is two thirds of the time required when GA is used and about one third of the settling time 
needed when PSO is used. Figure 5 and Figure 6 illustrates the advantages of using GSA to compensate P, 
and Q,, respectively. The settling time for the active power is at least 13 ms faster when GSA is considered. 
For the reactive power, GSA reduced the settling time by about 25% at least. 


V ar [P-0-] 


Vir [p-u.] 


Vie [p-u.] 


as i ] 10 10.5 11 
x 10 10.5 11 Time [sec] 


(a) (b) 


Figure 3. The rotor voltages when three-phase to ground fault is applied using GA, PSO and GSA: (a) v a 
and (b) iy 


GA 


10 10.5 11 
Time [sec] 


Figure 4. The rotor angular speed when three-phase to ground fault is applied using GA, PSO and GSA 
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Ps [p-u.] 


Py [p-u.] 


Ps [p-u.] 


10.5 
10 10.5 ll Time [sec] 


(a) (b) 


Figure 5. The stator active and reactive power when three-phase to ground fault is applied using GA, PSO 
and GSA: (a) P, and (b) Q, 


0.25 p 


Vor [p-u.] 
qr [Pe] 


> -005 


-0.15 F 


-0.25 L i | i i J 
10 11 12 13 14 15 16 


Time [sec] 


(a) (b) 


Figure 6. The rotor voltages when the reference speed is changed using GA, PSO and GSA: (a) var (b) Var 


4.2. Variable Reference Speed 
In this disturbance and as shown in Figure 6, the reference speed @m, ref Was initially set to 0.7 


p.u., after that, the speed is reduced to 0.35 p.u. in one second. Next, the reference speed stays constant for 
one second and finally increases linearly to the initial reference value. 


© [p.u.] 


Time [sec] 


Figure 7. The rotor speed with changing the reference speed using GA, PSO and GSA 
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From Table 2, a smaller ITAE is associated with the GSA as seen in Figure 6. The settling time of 
the output variables listed in Table 3 is damped faster by at least 0.38 s when GSA is used. The rotor voltages 
shown in Figure 7 clarifies that both GA and PSO give almost a similar response while GSA can reduce the 
maximum over shoot during the acceleration or the deceleration periods. 

The stator active and reactive power shown in Figure 8 clearly show the advantage of using GSA to 
damp the response much faster. Using GSA, the system requires 0.47 s to damp the reactive power while it 
takes 150% of this time when GA or PSO are adopted. 


-0.65 4 i i i ! i 
10 11 12 13 14 15 16 


10 11 12 13 14 15 16 Time [sec] 


(a) (b) 


Figure 8. The stator active and reactive power under variable reference speed using GA, PSO and GSA: 


(a) P, (b) Q, 


5. CONCLUSIONS 

In this paper, GSA is used to design optimal controllers for improving the transient response of a 
nonlinear DFIG-based WECS. The designed optimal controllers are tested under different common types of 
disturbances. The results are compared to the results obtained using the two well-recognized global 
optimization algorithms, which are GA and PSO, to show the effectiveness of using GSA to attain a global 
optimal solution of the design problem under study. A time-weighted absolute error (/TAE) is used in the 
objective function to stabilize the system and increase its damping when subjected to different types of 
disturbances. The comparative simulations confirm the fruitiness of the GSA and its effectiveness to damp 
the system faster than the GA and the PSO. Mathematically, the considered measured performance indices 
(ITAE and the settling time) show clearly the advantages of applying GSA over the other two algorithms for 
different disturbances. 
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